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Abstract. Morphological investigations and classification of Ge hut clusters forming 
the arrays of quantum dots on the Si(OOl) surface at low temperatures in the process 
of the ultrahigh vacuum molecular beam epitaxy have been carried out using in situ 
scanning tunnelling microscopy. Two main species of Ge hut clusters composing the 
arrays — pyramidal and wedge-shaped ones — have been found to have different atomic 
structures. The inference is made that shape transitions between pyramids and wedges 
are impossible. The nucleation probabilities of pyramids and wedges equal 1/2 at the 
initial stage of the array formation. The wedges become the dominating species as 
the amount of the deposited germanium is increased. A fraction and a density of the 
pyramids in the arrays are rapidly decreased with the growth of Ge coverage. 

The derivative types of the clusters — obelisks (or truncated wedges) and accreted 
wedges — have been revealed and investigated for the first time, they have been found 
to start dominating at high Ge coverages. The obelisks originate from the wedges as 
a result of their height limitation and further growth of trapezoid facets. The apexes 
of the obelisks are formed by sets of the parallel (001) ridges. 

The uniformity of the cluster arrays have been evidenced to be controlled by the 
length distribution of the wedge-like clusters. At low growth temperatures (360°C) 
nucleation of new clusters is observed during the array growth at all values of Ge 
coverage except for a particular point at which the arrays are more uniform than at 
higher or lower coverages. At higher temperatures (530°C) cluster nucleation has not 
been observed after the initial stage of the array formation. 

PACS numbers: 68.37.Ef, 81.07.Ta 
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1. Introduction 

1.1. Problem definition 

The development of processes of the controllable formation of germanium quantum dot 
(QD) arrays on the silicon surface as well as multilayer Ge/Si epitaxial heterostructures 
on their basis is a subject of significant and permanently increasing efforts for a number 
of years [H Ej [3l Hj primarily due to their potential applications in prospective devices of 
microelectronics and integrated microphotonics compatible with the monolithic silicon 
VLSI technology. Both high density of the germanium nanoclusters (> 10^^ cm^^) 
and high uniformity of the cluster shapes and sizes (dispersion < 10 %) in the 
arrays are required for many practically important application of such structures 

The main technique of formation of the germanium nanoclusters on the silicon 
surface is the molecular beam epitaxy (MBE) [2l [3]. A high density of the self- 
assembled clusters can be obtained in the MBE process of the Ge/Si(001) structure 
formation when depositing germanium on the silicon substrate heated to a moderate 
temperature (< 550°C)Iil In this case the lower the temperature of the silicon substrate 
in the process of the germanium deposition the higher the density of the clusters is 
at the permanent quantity of the deposited germanium [U |12l |13]. For example the 
density of the germanium clusters in the array reached 6 x 10^^ cm~^ at the substrate 
temperature during the deposition Tg^ = 360°C and the effective thickness of the 
deposited germanium layeil§| /ice = 8 A whereas the cluster density of only about 
2 X 10^^ cm~^ was obtained at Tgr = 530°C and the same value of Hq^ 

There are also different ways to increase the cluster density in the arrays. Thus, 
the authors of Ref. succeeded to reach the cluster density of about 9 x 10^^ cm~^ in 
the array using the pulsed irradiation of the substrate by a low-energy Ge"*^ ion beam 
during the MBE growth of the Ge/Si(001) heterostructure at Tgr = 570°C. 

Obtaining of the arrays of the densely packed Ge QDs on the Si(OOl) surface is an 
important task but the problem of formation of uniform arrays of the Ge clusters is even 
more challenging one. The process of Ge/Si(001) heterostructure formation with the 
Ge QD dense arrays and predetermined electrophysical and photoelectric parameters 
cannot be developed until both of these tasks are solved. The uniformity of the cluster 
sizes and shapes in the arrays determines not only the widths of the energy spectra 
of the charge carrier bound states in the QD arrays [5] but in a number of cases the 
optical and electrical properties of both the arrays themselves and the device structures 
produced on their basis [EJ. To find an approach to the improvement of the Ge QD 

X Remark also that lowering of the array formation temperature down to the values of < 450° C is 
required for the compatibility of the Ge/Si(001) heterostructure formation process with the CMOS 
device fabrication cycle [1] . It is another reason to decrease the temperature of all treatments starting 
from the Si surface preparation. 

§ I.e. the Ge coverage or more accurately the thickness of the Ge film measured by the graduated in 
advance film thickness monitor with the quartz sensor installed in the MBE chamber. 
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Figure 1. In situ STM micrographs of the Ge hut cluster arrays formed on the 
Si(OOl) surface in the process of molecular beam epitaxy at the substrate temperature 
Tgr — 360° C and different effective thicknesses of the deposited Ge layer (/iqc), the 
values of /iCe are (a) 6 A; (b) 8 A; (c) 10 A; (d) 14 A. Diagonals of the images are 
parallel to the (1 0) axes. 



array uniformity on the Si(OOl) surface it is necessary to carry out the morphological 
investigation of the clusters constituting the arrays and first of all classify them. 

At present, two genera of the self-assembled Ge clusters formed on the Si(OOl) 
surface are marked out — huts and domes. The former are smaller and faceted by 
{105} planes while the later are much larger and have more sophisticated faceting 
P, [131 [iHl [El [13 [18]. Our investigations of the densely packed Ge nanocluster arrays 
on the Si(OOl) surface [1] have shown that the composition of an ensemble of the hut 
clusters is by no means homogeneous — there are several species of the hut clusters 
different by their geometrical shapes as well as their behaviour in the process of the 
array formationjJJ] Examples of the Ge cluster arrays formed on the Si(OOl) surface at 
different Tgr and /ice are shown in Figs. [H and [2] which obviously demonstrate that the 

II Here we do not consider defects of the arrays. Other article will be devoted to their investigation. 
Now the available information about the morphology, structure and density of the array defects as well 
as their effect on the array parameters can be found in our reports [T] or [12] . 
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Figure 2. The same as in Fig. [U but for Tg,. = 530°C, Hqc- (a) 8 A; (b) 10 A; (c) 
11 A. 



arrays grown at low temperature always consist of a set of morphologically different hut 
clusters. Some of them have often been discussed in the literature since the classical 
letter hy Mo et al. [ISl [III [201 [HI 1221 [231 [21 [251 [26l [27] and some have not. 

This paper is devoted to the study of morphological differences of the Ge hut clusters 
formed in the process of molecular beam epitaxy at low substrate temperatures on the 
Si(OOl) surface and their classification derived from the revealed species differences. 

Writing about the classification of the hut clusters we have not got a goal to simply 
introduce a new terminology, as it might seem, although the latter is also proposedj^j 
The aim of classification is to sort out the hut clusters in accordance with their structural 
peculiarities which are much more important attributes than geometrical shapes. The 
difference of geometrical shapes does not necessarily imply the difference of atomic 

% Since the pioneering paper by Mo et al. [16j . a good few of descriptive and often confusing terms 
have been used in the hterature to designate two known types of the hut clusters. In this paper, we 
introduce a new strict stereometrical terminology to emphasize the structural difference between the 
cluster species and avoid muddle in the future. We shall name each species of the clusters in accordance 
with the denominations of the geometrical bodies which most accurately describe the shapes of the 
clusters. 
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structures. The latter may be identical. Hence, if so, it may be assumed that one type 
of clusters originates from another, as it is usually accepted "by default" (pastulated) 
for pyramids and "elongated" hut clusters in the literature [Ml [25| [26] . On the contrary, 
if the atomic structure of two species of clusters is different the shape transitions seem 
to be very unlike because of a probably high potential barrier to be overcome to change 
the cluster atomic configuration and symmetry. 

1.2. Formation of Ge hut clusters. A brief excursion into the historical domain 

It is well known [28j that in equilibrium Ge layers grow on the Si(OOl) surface following 
the Stranski-Krastanow mode [29l [30] . It means that initially Ge grows layer by layer 
until reaches a thickness of a few monolayers, then nucleation of the three-dimentional 
islands begins |2|. 

Because of lower surface energy, Ge when growing wets the Si surface (this enables 
tracing an etymology of the term "wetting layer"). Despite the mismatch of Si and 
Ge lattices, which is ~ 4.2 %, Ge atoms deposited on the Si(OOl) surface keep the 
correspondence with the Si atoms for several monolayers. The Ge layer surface is 
dimerized and (2x1) reconstructed [28]. Due to dimer buckling it comprises a mixture 
of the c(4 X 2) and p{2 x 2) structures [HI [32] which are manifested in STM images as 
characteristic antiphased and cophased zigzags. As Ge atoms arrive onto the surface, 
the compression along the dimer rows is relieved by arising dimer vacancies. Further 
ordering of the dimer vacancies in a nearly regular array of parallel trenches (or formation 
of the so called (2 x n) reconstruction [28j), which are perpendicular to the dimer rows, 
goes on to reduce the surface strain energy beyond the Ge coverage of 0.8 monolayer 
(thickness of 1 ML ^ 1.4 A) [33]. And at last, the formation of the dimer-row vacancies 
and a grating of the quasiperiodic (MxN) patched structure [321 [331 [311 [35] exhausts the 
ability of the dimer vacancies to accommodate the wetting layer to the strain increasing 
with the growth of quantity of the deposited Ge. When the Ge coverage exceeds 3 ML, 
the three-dimentional nanoislands (or Ge nanoclusters) start to nucleate on the surface. 
They are free of dislocations in the Si/Ge interface (coherent with the Si substrate), 
faceted and have base sides alined with two orthogonal (1 0) axes. At moderate growth 
temperatures the composition of the cluster arrays is bimodal: A part of clusters have 
shapes of regular square-based pyramids while others have rectangles in their bases 
[T6l [36] . Due to the shapes resembling huts both square-based and rectangular-based 
clusters are usually referred to as "hut" clusters. The hut clusters (coherent islands) 
were theoretically shown to be (under some conditions) more energetically stable than 
the strained films or dislocated islands [201 EZ] • Their appearance was also found to be 
kinetically favourable compared to the nucleation of dislocations [37j. 

The first STM observation of hut clusters was reported by Mo et al. in 1990 [H] 
(the term "hut" was introduced in that article). That letter presented an experimental 
investigation of a newly discovered metastable phase of Ge clusters which arose in the 
process of Stranski-Krastanow growth to relieve the increasing stress in the wetting layer 
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when Ge was deposited by MBE at moderate temperature on Si(OOl). As distinct from 
the macroscopic clusters, hut clusters were found to have {105} faceting on all sides. 
The first model of the cluster {105} faces according to which these facets consist of 
(001) terraces composed by pairs of dimers of the Ge(001)-(2 x 1) reconstructed surface 
was suggested in that communication (now this model is usually referred to as the PD 
model |36l [381 [39] ) . Huts were found to have "predominantly a prism shape (with canted 
ends), in some cases four-sided pyramids, with the same atomic structure on all four 
facets". The observed length of the huts was really huge (up to 1000 A) while their 
widths did not exceed 200 A (the aspect ratio often reached 10). The question was put 
for the first time in that article what caused the elongation and specific base orientation 
of huts. Unfortunately, no convincing answer to this question has been proposed thus 
far. 

It should be noted that since the letter by Mo et ai, pyramidal and "elongated" 
clusters have alway been considered in the literature as structurally identical ones 
different only by their base aspect ratio [211 [22l [23l [Ml [251 [26l [27] . The only argument 
for this assumption — the identity of faceting — does not seem to us to be very sohd. 
Stress relaxation via formation of the {105} faceted structures, such as islands or pits 
[2T] , appears to be energetically favourable and the structures are more stable compared 
to ones with different faceting [39] . This means that faceting itself cannot be considered 
as the only sign of belonging to some specific group of morphologically identical clusters. 
Total atomic structure of clusters is defined by both the structure of their facets and the 
configuration of their apexes. If the latter is different, the clusters should be regarded 
as members of different species. 

The sound counterargument to the assumption about the cluster identity, which 
is usually disregarded, was given, e.g., in Ref. [IQ]: "Elongated" clusters completely 
disappeared from arrays after annealing at 550°C for 600 s whereas pyramids and domes 
remained. Moreover, a commonly adopted pathway of the dome cluster formation is 
as follows: some "prepyramid" (we have never observed such formations in our low- 
temperature MBE experiments) pyramid dome Hlj . "Elongated" huts have 
never been met on this pathway. This led us to suggestion that "elongated" hutfl differ 
from pyramids not only by shapes but mainly by their atomic structure and, that is 
more important, by the genesis and the role in the array development. 

Nevertheless, the question about the mechanism of hut elongation was asked and 
required answering. Before long, the simplest and at first glance the most probable 
scenario was proposed and immediately adopted by the community. According to this 
scenario wedge arise by elongation of a pyramid due to the growth of one of its {105} 
facets [211 [221 [231 [211 [251 [261 [27] . This hypothesis would explain everything unless the 
necessity to explain why the symmetry of the pyramid is violated. Due to its shape the 
square-based regular pyramid seems to be stable enough, at least unless some exterior 
anisotropic agent affects it removing the degeneracy of its facets. Otherwise it is unclear 

Hereinafter we shall, as a rule, refer to them as wedges, wedge-like or wedge-shaped clusters. 
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why one facet gathers arriving Ge atoms to the prejudice of the rest three (or at least 
two) equivalent planes. 

The next often observed miracle is the formation of two closely neighbouring clusters 
separated by only a few nanometers, one of which is a pyramid while another is a wedge. 
This means that the acting agent violating the system symmetry is localized within a 
few nanometers around the growing wedge and does not affect the adjacent growing 
pyramid. At high enough coverages such clusters often start to coalesce — this in turn 
means that the agent previously violated the symmetry of one cluster of the pair has 
reached the growing pyramid but now does not affect it! 

It may be, however, that the symmetry is violated only once at the beginning of the 
pyramid elongation, then agent stops acting and the cluster grows on triangular facets 
faster than on trapezoid ones. It should be concluded now that, as distinct from the 
case of annealing [IQ], pyramids are much less stable than wedges in the process of Ge 
deposition, perhaps due to supersaturation by Ge atoms on the surface. This conclusion 
seems to be in agreement with experimental data, although the nature of the fluctuating 
agent stays unclear. 

Goldfarb et al. proposed that an asymmetry of the stress field due to the presence 
of {105} faceted pits may result in the cluster anisotropic lateral growth because of 
so-called equilibrium-driven elongation [26] . They found that in thick hydrogen rich Ge 
wetting layer (between 7.7 and 8.3 ML), which is formed in the process of gas-source- 
molecular-beam-epitaxy (GS-MBE)0 growth of Ge on Si(OOl) from GeH4, the strain 
at 690°C is reheved by formation of {105} faceted pits rather than islands [21], that 
is in agreement with the earlier conclusion by Tersoff and LeGoues [37] according to 
which pits always have a lower energy than islands of the same shape and equal size. 
(Goldfarb et al. observed also that at 620°C — and hence thinner wetting layer — only 
hut clusters arose [21].) Further, with the increase of Ge coverage, when the capability 
of pits to release the stress is exhausted hut clusters nucleate in the vicinity of pits 
[2T| [22| [25] . Acknowledging a model by Jesson et al., which explains instability of the 
hut cluster shapes (read "elongation") by nucleation and growth on the facets [27], as 
more common, the authors of Ref. [21] illustrate the elongation process by an example of 
cluster coalescence. We agree with Goldfarb et al. that such event sometimes happens 
and will present below a picturesque evidence of it. However, this explanation of the 
cluster elongation phenomenon of course by no means can be accepted as universal|il 

As mentioned above, investigations by Goldfarb with co-authors |26] eventually 
gave a weighty argument in support of the so-called equilibrium-driven elongation or, in 
other words, elongation governed by energy minimization. They considered the cluster 

t GS-MBE much more resembles chemical vapour deposition (CVD) than the solid source ultrahigh 
vacuum (UHV) MBE in which an atomic beam of Ge supplies the growing layer with Ge atoms rather 
than a flux of GeH4 or Ge2Hg. 

X Note, by the way, that Goldfarb with co-authors were first who reported the STM observation of 
incomplete trapezoid facets of the wedges in Ref. [22] . Until now, incomplete triangular facets have not 
been observed, however. 
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evolution model by Tersoff and Tromp developed for conditions when an isolated faceted 
strained Stranski-Krastanow island formed on a wetting layer grow in height slower than 
in lateral direction [20] . It may be shown in this case by minimizing the total energy per 
unit volume that there is a critical size of a pyramidal cluster (ao) and the cluster grows 
isotropically up to some point (eao), then anisotropic elongation in one direction starts. 
The authors of Ref. [26] indicate the discrepancy of the estimates made on the basis 
of the model by Tersoff and Tromp [20] with experimental observations. The model 
predicts eao ~ 100 nm (ao ~ 37 nm) whereas in opinion of the authors of Ref. [26] 
anisotropic elongation of huts starts at the nucleus size of ~ 5-8 nm [21j, if at all. 

Another model of the equilibrium driven growth was proposed by Li, Liu and 
Lagally for a two-dimentional rectangular island [261 US]- The model is as follows: 
Let dimensions of the island sides be s and /. Till its size is less than some critical 
value, the total energy is minimum for s = I. For a greater island its square shape 
becomes instable because of the strain: The island begins elongation in one of two 
degenerate orthogonal directions until the total energy reaches minimum at some value 
of arctan(s//) = 7r/4± A. To explain the anisotropic elongation the authors of Ref. [12] 
have to introduce anisotropy of edge energies to violate the symmetry of a square. 
According to this model island grows along the direction of the lower edge free energy 
in which both strain and edge energies are minimum. The advantage of this construct 
is that for strong anisotropy the elongation occurs at any size of the island. Goldfarb 
and co-authors extend this model to the case of a 3D faceted island with slowly growing 
height [26] . According to them when island elongates its perimeter grows faster than 
the area of the strained base in such a way resulting in more effective relaxation. 
Nevertheless, an origin of anisotropy stays an issue in this model. 

As mentioned above an explanation was however proposed in the same letter [2^. It 
was observed that huts interact with parent pits as well as with different adjacent ones, 
their lengthwise growth often starts at one pit and ends at the other. Sometimes the 
elongation is finished when a hut grows along the pit boundary and reaches its corner. 
The effect of the hut and pit interaction was analyzed by finite element technique and 
found to cause hut elongation. The pit dimensions used in calculations were 10 x 10 nm 
and the wetting layer was as thick as 2nm. The inference was made in the article that 
the following hut evolution scheme takes place: as soon as a stable critical nucleus is 
formed [25] it grows in the energetically favorable direction along its mutual boundary 
with the parent pit until reaches the pit corner or attains an equilibrium s/l ratio. This 
finishes the first phase of cluster elongation. The second phase implies that the cluster 
continues the equilibrium growth in the perpendicular direction (from the pit) until 
either an equilibrium s/l ratio is established or impingement to other pit happens. 

These observations and reasoning maybe correct for the particular case of GS-MBE 
were propagated to all deposition methods. We have solid objections to it. First of all, 
only relatively rarefied arrays on very thick hydrogen-rich wetting layers were grown by 
GS-MBE and investigated in Ref. [26]. Dense arrays obtained by UHV MBE grow in 
absolutely different manner (Figs. [1] and [2]). Secondly, as it is seen from Figs. [1] and [21 no 
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pits are available in the wetting layer grown by UHV MBE at moderate temperatures. 
Thirdly, according to our observations, which are as a rule carried out with atomic 
resolution, neither pits nor steps are required for cluster nucleation and growth by UHV 
MBE (Figs. [Hand [21 Ref. |43]). And at last, the described tricky evolutions of hut would 
leave their imprints on the cluster shape and structure which would be seen in a good 
microscope. Unfortunately, this is not the case and no marks of the above evolutions 
are seen even at atomic resolution. 

A model competing with those of equilibrium-driven elongation is usually referred 
to as kinetically driven elongation [23l [2ll [H] . In this model, one of the pyramid facets 
begins to grow by chance due to a fluctuation. Coverage by a monolayer decreases 
the likelihood of the next attachment of a Ge adatom to newly appeared trapezoid 
facets whereas the probability of attachment to triangular ones remains constant. Due 
to increase of areas of the trapezoid facets and consequently the barrier of adatom 
attachment, further elongation in the randomly chosen (10 0) direction goes on with 
decreasing probability of Ge atom attachment to the trapezoid facets and hence with 
decreasing rate of the in-height growth in comparison with the rate of the longitudinal 
one. Unfortunately this model appeared to disagree with experiments. Observations 
of Goldfarb et al. |26j indicate that in case of interaction of pits and huts the two 
(10 0) directions of elongation are not equally probable as it follows from the model of 
kinetically driven elongation. Our elaborated observations, which are presented below, 
do not support this model either. 

Summarizing this brief historical review we would like to emphasize the following. 
To date, main milestones on the pathway of the Stranski-Krastanow growth of Ge film on 
Si(OOl) at moderate temperatures (< 550°C) from the pure silicon surface to hut clusters 
appearance are recognized as consecutive steps of the strain relief. They are following: 
(2 X 1) ^ (2 X ^ (M X N) —>■ huts. However, despite the efforts made there 
is no clarity in the issues of hut nucleation and further transformation. The process 
of hut nucleation appeals for detailed experimental investigation at different growth 
conditions by instruments assuring atomic resolution and using different deposition 
methods. Formation and longitudinal growth of huts have not been understood thus 
far. Elongation of pyramids has never been unambiguously observed in experiments. 
The best of the two theoretical models explaining the elongation of wedges also has 
not been chosen yet. Issues of evolution of the cluster arrays during Ge deposition has 
been passed over by researchers too. There has been no systematic investigations of the 
stages of this very important and complicated process presented in literature. A final 
phase of the array growth cycle — growth at high coverages and transition to the 2D 
mode — has never been in focus of investigations. So, we can conclude now that despite 
the widely adopted standpoint, investigations on the discussed problem are still far from 
completion. 

The above analysis of the literature puts a number of new questions the most 
obvious of which are as follows. Firstly, it is unclear whether the structure of pyramids 
and wedges is identical and whether they belong to the same morphologically uniform 
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class. Secondly, does the structure of pyramids and wedges coincide at the moment of 
nucleation or have they got different structures aheady at the stage of occurrence and 
consequently different nuclei? Thirdly, are their roles the same in the array formation 
and evolution? Fourthly, what are the driving forces of their evolution during the array 
growth? Why the results of the evolution are different for pyramids and wedges? And 
at last, is it possible to control the cluster evolution to obtain uniform and defectless 
dense arrays suitable for the industrial applications? 

In the current paper, we shall study a part of the above issues focusing mainly on 
the morphology of hut clusters and partly on the growth cycle of the dense arrays. 

2. Experimental 

The experiments were made using an integrated ultrahigh vacuum system based on the 
Riber EVA 32 molecular beam epitaxy chamber coupled through a transfer line with 
the STM GPI-300 |15] scanning tunnelling microscope (STM). This instrument enables 
the STM study of samples with atomic resolution at any stage of Si surface cleaning 
and MBE growth. The samples can be consecutively moved into the STM chamber 
for analysis and back into the MBE chamber for further processing never leaving the 
UHV ambient and their surfaces stay atomically clean over whole the experiment. A 
procedure of the sample preparation was as follows: initial substrates were 8x8 mm^ 
squares cut from B-doped CZ Si(lOO) wafers (p-type, p = 12 ficm). After washing 
and chemical treatment following a standard procedure described elsewhere (see e.g. 
Refs. [H HE]) the silicon substrates mounted on the molybdenum STM holder and 
clamped with the tantalum fasteners were loaded into the airlock and transferred to the 
preliminary annealing chamber where outgassed at the temperature of around 565° C and 
the pressure of about 5 x 10~^ Torr for about 24 hours. After that the substrates were 
moved for final treatment into the MBE chamber evacuated down to about 10^^^ Torr. 
There were two stages of annealing in the process of substrate heating in the MBE 
chamber — at ~ 600°C for ~ 5 min. and at ~ 800°C for ~ 3 min. (Fig. [3]). The final 
annealing at the temperature greater than 900°C was carried out for nearly 2.5 min. 
with the maximum temperature of about 925°C (~ 1.5 min.). Then the temperature 
was rapidly lowered to about 750° C. The rate of the further cooling down was around 
0.4°C/s. The pressure in the MBE chamber enhanced to nearly 2 x 10^^ Torr during 
the process. 

The surfaces of the silicon substrates were completely purified of the oxide film as 
a result of this treatment. The high-order (8 x n) surface reconstruction described in 
Ref. ^\ was always revealed by the STM on the deoxidized substrates whereas the 
reflected high energy electron diffraction (RHEED)[§| patterns obtained from the same 
purified surfaces always corresponded to either (2 x 1) or (4 x 4) surface reconstruction 
|48] . This observation is in a good agreement with the model brought forward by 

§ RHEED is usually applied in the MBE vessels for monitoring the surface perfection [3] e.g. during 
the deoxidizing process. 
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Figure 3. A diagram of tlie final tliermal treatment of tlie Si substrates during surface 
deoxidizing in the MBE chamber. 



US in Ref. |17] as well as with the generally accepted opinion based on the RHEED 
measurements that the (2 x 1) reconstruction is formed on the Si(OOl) surface due to 
deoxidization in an MBE chamber [5j. 

Germanium was deposited directly on the purified silicon surface from the source 
with the electron beam evaporation. The rate of Ge deposition was about 0.15 A/s, 
/iQe was varied from 6 A to 14 A for different samples. The deposition rate and the 
effective Ge film thickness Hq^ were measured by the graduated in advance XTC film 
thickness monitor with the quartz sensor installed in the MBE chamber. The substrate 
temperature Tgj. was 360° C or 530° C during the process. The pressure in the MBE 
chamber did not exceed 10^^ Torr during Ge deposition. The rate of the sample cooling 
down to the room temperature was approximately 0.4°C/s after the deposition. 

The samples were heated by Ta radiators from the rear side in both preliminary 
annealing and MBE chambers. The temperature was monitored with chromel-allimel 
and tungsten-rhenium thermocouples in the preliminary annealing and MBE chambers, 
respectively. The thermocouples were mounted in vacuum near the rear side of the 
samples and in situ graduated beforehand against the IMPAC IS 12-Si pyrometer which 
measured the sample temperature through chamber windows with an accuracy of 
±(0.3% r°c+ 1°C). 

The atmosphere composition in the MBE camber was monitored using the 
SRS RGA-200 residual gas analyzer before and during the process. 

After Ge deposition and cooling, the prepared samples were moved for analysis into 
the STM chamber in which the pressure did not exceed 10~^° Torr. The STM tip was ex 
situ made of the tungsten wire and cleaned by ion bombardment [49j in a special UHV 
chamber connected to the STM chamber. The images were obtained in the constant 
tunneling current mode at the room temperature. The STM tip was zero-biased while 
the sample was positively or negatively biased for empty or filled states mapping. 

The WSxM software [50] was used for processing of the STM images. 
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Figure 4. STM micrographs of tlie pyramidal Ge clusters; completely shaped clusters: 
Tg, = 530°C, /iGc = 11 A (a); T^, = 360°C, /iqc = 10 A (b); a cluster with unfinished 
facets (marked by an arrow): Tgr = 360°C, /iqc = 14 A (c); a structure of vertex, sides 
and edges (d) and a structure of a nucleus of 1 monolayer high over the wetting layer 
(e): Tgi- = 360°C, /igc = 6 A (corresponding features on toppings are marked by the 
same numerals). 
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3. Classification 

3.1. Main species of the clusters 

3.1.1. Pyramidal and wedge-shaped clusters 

As mention above an array of tlie self-assembled germanium hut clusters formed on 
the Si(OOl) surface consists of a set of morphologically different clusters. All the clusters 
have the edges of bases oriented along the (100) directions in common. Yet in spite of 
the apparent variety of the cluster forms (Figs. [1] and W) an analysis of the STM images 
gives evidence that only two main species of the hut clusters exist — ones having square 
bases and shapes of the regular pyramids and those with rectangular bases which have 
shapes of the wedges. We have already cited the paper by Mo et al. [H], in which 
both species of the hut clusters were described for the first time, as well as a number of 
publications which investigated the details of their formation [201 HH [22], [231 [211 [251126] . 
Unfortunately their structure has not been yet clearly visualised and identified. 




X: S.lnm (c) 

Figure 5. STM images of the wedge-like Ge clusters; a cluster with an unfinished side: 
Tgr = 530° C, /iGo = If A (a); the entirely formed cluster: Tgr = 360° C, /igc = f A (b); 
a structure of the ridges (two closely neighbouring clusters), the shifted features are 
marked by the arrows: Tgr = 360°C, /iGo = 8 A (c). 



Let us dwell on the description of each species of the clusters in more detail. 
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Typical high resolution STM images of the pyramidal clusters in the arrays with 
different growth parameters are shown in Fig. HI A regular shape of the clusters is 
clearly seen in the pictures (a) and (b) presenting the arrays obtained at Tgr = 530°C, 

= 11 A and Tg^ = 360° C, Kq^ = 10 A, a fine structure of the faces is resolved as 
well as the (M x N) patched structure of the wetting layer [32], [33], [311 EH] . Lines of 
the solidifying steps are revealed for the first time on the cluster incomplete faces in the 
image (c) of the array obtained at T^, = 360°C, /ice = 14 A (one of them is marked 
by arrow in the image). Fig. [D^d) shows a vertical view of the small pyramid grown at 
Tgr = 360°C, /iGe = 6 A and having only 5 monolayers height over the wetting layer. 
A fine structure of the pyramid vertex and edges as well as the stepped structure of its 
{105} facets are resolved in detail in the image (d). And at last, the same structure as 
that seen in Fig. [4](d) on the pyramid vertex is clearly resolved in the image (e) of the 
pyramid nucleus (1 monolater high over the wetting layer) situated on the block of the 
Ge (M X N) surface jj 

Fig. [5] demonstrates STM images of the wedge-like clusters. Being the hut clusters 
they are bounded by {105} planes, i.e. their heights are to the base widths as 1:10. 
A distinctive feature of this species of the clusters is that their base lengths are not 
connected with cluster heights and are rather random. To some extent, the base lengths 
of the wedge-like cluster depends on its nearest neighbours. Nevertheless, it is impossible 
to confidently point out the factors which affect the lengths of the Ge wedges based upon 
the available data. Reasoning from the results of the STM image analysis it may only 
be asserted that their base length-to-width ratio is distributed randomly and rather 
uniformly in the interval from a little greater than 1 to more than 10. 



5 



4 



3 

E 
a 

N 

2 








10 20 30 40 50 60 
X, nm 

Figure 6. Profiles of tlie neigtibouring wedge (the left one, taken along the short side 
of base) and pyramid (the right one) shown in Fig. [U^a) . 

II Note also that similar configuration of the pyramid apex can be discerned by an attentive observer 
in Fig. [H^b) which presents an image of a "ripe" pyramid in the well developed array. 
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As long as the factors determining the base lengths of the wedges are unclear 
it is unknown if it is possible to control the array growth parameters in such a way 
to minimize the dispersion of the base lengths of the wedge-like clusters. Moreover, 
the factors governing the nucleation of either pyramidal or wedge-like clusters of very 
different lengths on the wetting layer are also vague. The words about the role of 
the strain in the wetting layer themselves explain nothing. It stays unclear why the 
clusters with different shapes (or even symmetry) may arise in a very close vicinity 
to one another. Perhaps the nucleation process and consequently the strain field (its 
distribution and local symmetry) are controlled by the underlying Si surface of the Si/ Ge 
interface — its reconstruction before and defect (in particular, vacancy) distribution 
before and in the process of Ge deposition — especially in the case of the low-temperature 
processes. This is one of the reasons to consider the surface pre-growth treatment as a 
key to the controllable Ge/Si heterostructure formation process. That is why the in situ 
studies of the Si surface on the atomic scale immediately before the Ge deposition as 
well as the investigations of its influence on the deposited Ge layer should be considered 
as a task of high importance. 

Fig. [5](c) demonstrates the fine structure of ridges of two close wedge-like clusters 
(Tgr = 360°C, /iGe = 8 A). It is interesting that one can descry the same configuration 
of the ridge in the first published image of the hut cluster (see Fig. 2 in Ref. p6]). 

The STM images of fine structure of the vertexes and the ridges, similar to those 
shown in Figs. IH^d, e) and[5]^c), helped us propose structural diagrams of both species of 
the clusters [l3l [51]. It is clearly seen from the images that the fine structures of apexes 
of the clusters are different. The features in the uppermost parallel rows on the ridges 
of the wedge-like clusters are shifted with respect to one another. (They are marked 
by the rows of the shifted arrows in the STM image.) We interpret these features as 
Ge dimer pairs in accordance with the simple structural model of the hut cluster facets 
(PD model) proposed by Mo et al. [16j. Similar features in the images of vertexes of the 
pyramids are gathered in the straight rows. The difference of symmetry of the pyramid 
vertex from the symmetry of the elementary unit of the wedge ridge is distinctly evident 
when compared Figs. |l](e) and[5](c). 

As we have already mentioned above, the difference of the atomic structure causes 
a ban of shape transitions between the pyramidal and wedge-shaped hut clusters which 
were intensively discussed in the literature [211 [25], [26]. In addition, particular nuclei 
should be sought for pyramidal and wedge-shaped clusters. The question appears as 
well, why two structurally different species of hut clusters arise on the wetting layer. 

Fig. [6] shows the cross section profiles of the adjacent wedge-like and pyramidal 
clusters presented in Fig. [D^a). Both clusters are seen to have equal ratio of the base 
width to the cluster height close to 10. The base sides of the pyramidal cluster and 
the base length of the wedge-like one are nearly equal whereas the base width of the 
wedge-like cluster is by about 1.6 times less than the sides of the pyramid base. It is 
a common rule which is not affected by the length of a particular wedge-like cluster: 
The pyramidal clusters are usually higher than the wedge-like ones if /iQe is high enough 
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(see also, e.g., images (a) to (c) in Fig. H]). This observation may be explained if height 
limitation is assumed to occur for the wedge-like clusters at Hq^ greater than some value 
and not to occur for the pyramids. We do observed such limitation for the wedge-shaped 
clusters (see below) and did not succeed to fix the height limitation of the pyramidal 
onesJ 
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Figure 7. The density of tlie Ge clusters in the arrays (a) formed at Tgr — 530° C 
(designations: □ marks the pyramids, A corresponds to the wedges, O means the 
total value) and Tgr — 360°C (■ designates the pyramids, ▲ means the wedges, 
• is the total density). The fraction of the pyramidal and wedge- like Ge clusters 
in the arrays (b), Tgr = 530°C ( □ marks the pyramids, A designates the wedges), 
Tgi- — 360°C (■ corresponds to the pyramids, ▲ designates the wedges). Both graphs 
are plotted vs /iqc- 



3.1.2. Density and fraction of each species of clusters in the arrays 

It was observed that the wedge-like and pyramidal clusters are different not only 
in their atomic structure and geometrical shapes. The wedge-like clusters dominate in 
the arrays formed at low temperatures, and their fraction grows with the growth of /iQe 

(Fig. CD. 

Fig. [Tl^a) plots the dependence of the cluster density on Hqc for different clusters 
in the arrays grown at 360°C and 530°C. It is seen that for Tgr = 360°C the density 
of wedges rises starting from 1, 8 x 10^^ cm~^ at the beginning of the three- 

dimensional growth of Ge (the estimate is obtained by data extrapolation to /ice = 5 A) 
and reaches the maximum of ~ 5 x 10^^ cm~^ at /iQe ~ 8 A, the total density of 
clusters at this point D^, ~ 6 x 10^^ cm~^ is also maximum. Then both and 
Dy, slowly go down until the two-dimensional growth of Ge starts at Hq^ ~ 14 A and 

^ Perhaps, there is no height limitation for pyramids and namely they give rise to large clusters affecting 
the properties of the Ge/Si(001) heterostructures and classified by us as one of the types of defects of 
arrays (see report [I] or an article on defects of arrays [19^). 
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Dy. ~ -Dw ~ 2 X 10^^ cm~^ (the contribution of pyramids -Dp to -Ds becomes negligible — 
~ 3 X 10^° cm~^ — at this value of /ice)- The pyramid density exponentially drops for 
Tgr = 360°C as the value of /ice grows (-Dp ^ 5 x 10^^exp{— 2.0 x lO^/ice}, ^Ge is 
measured in centimeters). The maximum value of -Dp 1.8 x 10^^ cm^^ obtained from 
extrapolation to /ioe = 5 A coincides with the estimated initial value of -D^. 

For Tgr = 530° C, the total density of clusters exhibits the same trend as for 
Tgr = 360°C, D^^7 X W^exp{-A.3 x W Hq^ [cm]}. The maximum (initial) value of 
-Ds is estimated as 8 x 10^^ cm^^ by extrapolation to Hqc = 5 A. 

The graphs of cluster fractions in the arrays versus /iqc are presented in Fig. El^b). 
For Tgr = 360°C, portions of pyramids and wedges initially very close (~ 50% at 
^Ge ~ 5 A) rapidly become different as Hq^ rises. The content of pyramids monotonically 
falls. The fraction of the wedge-like clusters is approximately 57% at the early stage of 
the array growth {hce = 6 A) and becomes 82% at /ice = 8 A. At further growth of 
the array, the content of the wedges reaches the saturation at the level of approximately 
88% at hce = 10 A. 

At moderate values of Hq^, the proportion of pyramids to wedges for Tgr = 530° C 
was found to be nearly the same as for Tgr = 360°C. The content of the pyramidal 
clusters in the array is about 20% at /iqc = 8 and 10 A. 

The inference may be made from this observation that contrary to the intuitively 
expected from the consideration of symmetry, the wedge-like shape of the clusters is 
energetically more advantageous than the pyramidal one, and the more advantageous the 
more Ge atoms (and the more the number of atomic layers) constitute the cluster. The 
probability of nucleation appears to be close to 1/2 for both wedges-hke and pyramidal 
clusters at the initial stage of the array formation and low growth temperatures. Then, as 
the array grows, the formation of pjnramids becomes hardly probable and most of them, 
which have already been formed, vanish whereas the nucleation and further growth of 
wedges continues (Fig. [1]). The Ge pyramides on the Si(OOl) surface appear to be less 
stable species than the wedges and in accordance with the "bourgeois principle" ("the 
survival of the fittest" ) they loose their substance in favour of the wedge-like clusters. 

At higher temperatures, no nucleation of new clusters was observed in the process 
of the array growth (Fig. [2]). The "bourgeois principle" decreases the cluster density in 
the arrays and increases their sizes despite the species they belong to. 

It should be noted that the above analysis demonstrates that pyramidal and wedge- 
like clusters are really different objects which belong to different cluster species rather 
than the varieties of the same structurally uniform species as it is usually postulated in 
the literature [lEl EH El EHl [26] . 

Remark also that at Tgr = 360°C and the flux of Ge atoms d/ice/dt = 0.15 A/s, 
the point /ice = 10 A is particular. Not only the fraction of pyramids saturates at this 
point but the array in whole has the most uniform sizes of the clusters composing it 
(Fig. [1]). This is concluded by us not only on the basis of analysis of the STM images of 
the Ge/Si(001) arrays but also from the data of the Raman light scattering by the Ge/Si 
heterostructures with different low-temperature arrays of Ge quantum dots [52l [53] . We 
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refer to such arrays as optimal. 

A qualitative model accounting for the presence of the particular point at the low- 
temperature array growth is simple. The case is that at low enough temperatures of 
the array growth, the new Ge cluster nucleation competes with the process of growth 
of earlier formed clusters. The height of the clusters (at least the dominating wedge- 
like ones) is observed to be limited by some value governed by Tg^. At small /ice, Ge 
clusters are small enough and the distances between them are large enough compared 
to the Ge atom (or dimer) diffusion (migration) length on the surface for nucleation of 
new clusters on the Ge wetting layer in the space between the clusters (Fig. [TJ^a, b)). 
At Hqc = 10 a and the above d/iGc/dt value, the equilibrium of parameters (cluster 
sizes and distances between them, diffusion length at given temperature, Ge deposition 
rate, etc.) sets in, the rate of new cluster nucleation is decreased and the abundant Ge 
atoms are mainly spent to the growth of the available clusters (Fig. [11(c)). After the 
clusters reach their height limit and in spite of it, Ge atoms continue to form up their 
facets. As soon as most of the clusters reach the height limit, nucleation of new clusters 
becomes energetically advantageous again and the nucleation rate rises. A second phase 
of clusters appears on the wetting layer and fills whole its free surface as /ice is increased 
(Fig. [D^d)). Further increase of Hq^ results in two-dimensional growth mode. It is clear 
now why the array is the most homogeneous (optimal) at Tgr = 360° C and /iQe = 10 A 
whereas the dispersion of the cluster sizes is increased at higher and lower values of Hq^ 
because of the small clusters included in the array. It is clear also that the optimal array 
will appear at different value of /iQe when Tgr or d/iQe/dt are different. 




Figure 8. STM micrographs of the Ge wedges with two ridges (obehsks), a general 
view of the cluster (a) and a top view of the ridges (b); Tgr = 360°C, /iqc = 14 A. 

In addition, some threshold value of Tgj. must exist beyond which the cluster growth 
process always dominates and the nucleation of clusters happened once will never be 
repeated. An example of such arrays formed at Tgr = 530° exceeding the threshold value 
is given in Fig. O 
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X, Y, nm 

Figure 9. Profiles of the Ge obelisk presented in Fig. [5Ja) measured across the ridges 
along the base width (the solid line) and over the left ridge along the base length (the 
dashed line). Here, as distinct from Fig. [SJ X and Y mean respectively axes directed 
along the cluster base width or along its length. Height is as usually counted from the 
wetting layer level (Z = 0). 

3.2. Derivative species of clusters 

3.2.1. Obelisks (truncated wedges with two ridges) 

Except for the described above main species of Ge hut clusters, different clusters 
are also formed on the Si(OOl) surface which cannot be classified as independent species 
because they originate from the wedge-like clusters but have specific shapes, particular 
formation mechanisms, maybe peculiar properties and hence should be marked out in 
a separate but derivative species. Fig. [8] shows clusters related to one of the species 
of derivative clusters — truncated wedge-like clusters with two ridges or obelisk-shaped 
clusters {T^ = 360°C, hce = 14 A). 

Profiles of the Ge obelisk shown in Fig. [Uj^a) taken along the short and long base 
sides are presented in Fig. [9l Although these clusters are the huts and have a slope of 
the facets ~ 11.3° the ratio of the cluster height to its base width is ~ 0.06. 

Fig. [lO] presents the STM images of Ge wedges with two ridges explaining the way 
of their formation. The double ridges are seen to arise as a result of building of the 
trapezoid facets in the process of the cluster growth. Initially the clusters had the 
only ridge being ordinary wedge-like ones. Wedges with double ridges arise because a 
limit value of the cluster height exists for any temperature Tgi.. If the limit is reached 
the further growth of the cluster always goes on by building its trapezoid facets and 
increasing its width. 

This species of clusters dominates in the arrays at high values of /iQe which depend 
on the value of Tgr (see Figs. [Hand [2]). 

A unique illustration of the process of the trapezoid facet growth is shown in 
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X: 32.8mn (c) X: 80.3nm (d) 




Figure 10. STM images of the Ge obelisks illustrating the mechanism of their 
formation; Tg,. = 360°C, /iQe = 14 A (a, b, c); Tg, = 530°C, /ice = 8 A (d) and 
10 A (e, f). 

Fig.[TO](c). Several (from four to six) incomplete (001) terraces are seen near the bottoms 
of the clusters (the dimer pairs are distinctly resolved, the arrows show these new growing 
facets in the STM image). These incomplete faces are seen to repeat the shapes of the 
former faces on which they grow and which are also incomplete. The highest terraces 
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start to grow before the lower ones finish the completion process, as well as new facets 
nucleate before the old ones finish the growth. It is also observed that the new facets 
nucleate far from the corners adjoining the base sides. As a result, a compound structure 
of the trapezoid facets is formed which can be seen in the truncated wedges represented, 
e.g., in Fig. 111(c). We propose the reader to compare the above description with the 
speculations by Jesson et al. which draw a different picture of the facet growth |23j . 
First of all, they consider a triangular face as a preferential site of a new facet nucleation 
explaining in such a way elongation of clusters. Then, according to their model, the 
facets nucleate in the corners rather than somewhere else, etc. 

We would like to remark that we have not succeeded to observe the growth of the 
triangular faces at the temperatures as low as 360°C. Nevertheless we observed this 
process at Tgr = 530°C. Figs. fTOT d-f) demonstrate this phenomenon. The peculiarities 
of the process are as follows: A new facet formation takes place when the cluster has 
already reached its hmit height and its additional trapezoid facets are well developed. 
The growing triangular facets are clearly observed on only one side of the clusters. The 
triangular faces can nucleate both far from the bottom corner (Figs. [TOl(d)) and close 
to the corner (Figs. fTOVe. f)). The growing (incomplete) faces replicate the shape of 
initial facet even if the latter is complex (composed by intersecting triangles due to 
the developed additional trapezoid faces, see Figs. [TOlfe. f) in which the growing faces 
on the short sides of the truncated wedges are shaped by two combined triangles). 
It can be deduced from these observations that the described process of formation of 
new triangular facets is different from that resulting in the discussed above significant 
elongation of the wedge-shaped clusters at earlier stages of the cluster growth. 

Now we would like to attract the reader's attention to the observed in Figs.[TOl(c. e, f) 
formation of the so called "square based clusters" from the wedge-like ones which is 
caused by the extensive growth on the trapezoid facets and the cluster widening. In 
Fig. fTOTc). the nearly "square based clusters" are seen to be formed because of the 
widening of the wedges with two ridges. They resemble truncated pyramids but actually 
preserve the structure of the wedge. The nearly "square based clusters" are also seen in 
the upper left corner of Fig. [TOT e) and upper right corner of Fig. [TOT f). These clusters are 
formed of the formerly wedge-like ones by successive addition of new incomplete facets. 
Their faces are complex and their shapes are far from the shape of an ideal regular 
pyramid. Certainly, their structure stays that of the wedge. The genuine pyramidal 
cluster revealed in the upper right corner of Fig. fTOl d) grows nearly uniformly on all 
four its triangular faces (compare also with Fig.Hl^c)). We would like also to indicate the 
formation of serial incomplete faces resolved on the sides of both pyramidal and wedge- 
shaped clusters presented in Fig. [TOT d) (even the dimer pairs on the "parallel steps" are 
seen). This process do transform the shape of the clusters and may create nearly "square 
based clusters" from the wedges as well as so called "rectangular based clusters" from 
the pyramids (the latter process may be fancied, e.g., if a pyramid is closely surrounded 
by its neighbours from all sides except for one and have some room for elongation only 
in one direction). Of course, such transformed clusters are always "truncated", have 
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a complex "stepped" structure of successive incomplete facets and atomic structure of 
apexes characteristic for their precursors, as it is seen in the presented STM images. 

3.2.2. Accreted wedges 

In Fig. [m the STM images of accreted together wedge-like Ge clusters are shown — 
accreted clusters which formed ones with two ridges (Tgr = 360°C, Hq^ = 8 A), accreted 
clusters which gave rise a F-like one (Tgr = 360°C, /ice = 10 A), clusters accreted 
approximately at half width formed an extended one with zigzag on the ridge and the 
facets (Tgr = 360°C, Hqc = 10 A) are depicted. 




Figure 11. STM images of the accreted together Ge wedges; the accreted clusters 
borne the clusters with two ridges, Tgr = 360°C, /iqc = 8 A (a); the accreted clusters 
given rise to the F-shaped cluster, Tgr = 360°C, /igc = 10 A (b); the merged at about 
half-length wedges formed an extended wedge-like cluster with zigzag on the ridge and 
the sides, Tgr = 360°C, hce = 10 A (c). 

Like the obelisks the accreted clusters cannot be classified as independent species 
because they also originate from the wedge-like clusters. However, they also have 
specific shapes and probably peculiar properties and consequently like obelisks should 
be separated in a special but derivative species. It was found from the analysis of the 
STM images obtained at different stages of the array formation that the nuclei of such 
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clusters are situated at the distance of < 32 A from one another at the initial stage of 
the array formation in the applied array growth conditions. 

From physical viewpoint, these clusters are obviously single and integral objects, 
the properties of which may be different from the properties of the usual wedge-like 
clusters. Their influence upon the properties of the arrays in whole is in prospect of 
further investigations. 

The coalescence of clusters is also seen in Fig. Wc)- The truncated wedges are 
merged or even completely absorbed by the growing pyramid (the pyramids are always 
greater than the wedges). Such process usually takes place at high values of Hq^ just 
before the beginning of the two-dimensional growth. 

4. Conclusion 

Summarising the above we would like to emphasise the central ideas of the paper. 

Morphological investigations and classification of Ge hut clusters forming the arrays 
of quantum dots on the Si(OOl) surface at low temperatures in the process of the 
ultrahigh vacuum molecular beam epitaxy have been carried out using in situ scanning 
tunnelling microscopy. The study reported in the paper was made in view of the 
necessity to controllably produce highly uniform and very dense arrays of Ge quantum 
dots at low temperatures in the process compatible with the CMOS one. Although this 
task is still far from the solution an important step is made in understanding the object 
properties to be controlled. 

The Ge/Si(001) system appeared to be much more sophisticated than it seemed to 
most of the researchers, and the knowledge about it which is present in the literature 
now seems to be very deficient and sometimes incorrect. This seems to be the main 
cause of failure for the last two decades to develop electronic or photonic devices on the 
basis of ensembles of Ge quantum dots on the Si(OOl) surface. 

Analysis of the high quality STM images which can be obtained only using an 
integrated high resolution UHV STM-MBE instrument allowed us to introduce a 
new classification of germanium hut clusters formed on the Si(OOl) surface. The 
hut clusters were found to be subdivided into four species, two of which are basic 
and structurally different — the wedge-like and pyramidal clusters — and the rest are 
derivative — the obelisk-shaped and accreted wedge-shaped clusters. The conclusion 
was made that shape transitions between pyramids and wedges are prohibited. The 
nucleation likelihoods of pyramids and wedges appeared to equal 1/2 at the initial 
stage of the array formation. The wedge-like clusters were observed to quickly become 
the dominating species in the arrays while the pyramidal clusters were found to 
exponentially rapidly disappear as the arrays grow. 

The derivative types of the clusters have been found to start dominating at high Ge 
coverages. The obelisks originate from the wedg result of their height limitation 

and further growth of trapezoid facets. The apexes of the obelisks are formed by sets 
of the parallel (001) ridges. 
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At low growth temperatures (360°C) nucleation of new clusters is observed during 
the array growth at all values of Ge coverage except for a particular point at which 
the arrays are more uniform than at higher or lower coverages. At higher temperatures 
(530°C) cluster nucleation has not been observed after the initial stage of the array 
formation. 

The growing trapezoid and triangular cluster facets were visualised. The 
peculiarities of the facet completion were described. It was shown that the growth 
of incomplete facets results in a complex structure of the growing hut clusters. 

It was shown also that the uniformity of arrays is governed by the lengths of the 
wedge-like clusters. This parameter is hardly controllable as distinct to the cluster 
width which is bound to the cluster height and hence is much more predictable. The 
cluster lengths are now absolutely unpredictable. Moreover the origin which determines 
their values is unknown at present. This difficulty requires extensive investigations 
and intensive efforts to be overcome. But it worth while doing because both stochastic 
(disordered) and the most promising artificially ordered or self-arranging dense arrays of 
self-assembled Ge clusters on the Si(OOl) surface [51] are equally required to be uniform 
and equally subjected to effect of the Ge wedge length unpredictabihty. 
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